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High-frequency  (HF)  radio  signals  propagated  to  large  distances  by 
ionospheric  refraction  remain  of  practical  interest  as  a  means  for  commu¬ 
nication  and  remote  probing  (e.g.,  over-the-horizon  radar).  Scatter  from 
ionospheric  irregularities  can  strongly  affect  such  signals.  An  under¬ 
standing  of  the  effects  of  such  scatter  therefore  is  necessary  if  the 
performance  of  systems  that  make  use  of  such  signals  is  to  be  determined. 

A  useful  characterization  of  scatter  effects  at  HF  is  available  for 
the  natural  environment,  but  it  has  been  attained  almost  entirely  empiri¬ 
cally.  Our  knowledge  of  the  propagation  environment  following  high-altitude 
nuclear  detonations  is  more  fragmentary.  High-frequency  ionosonde  data 

taken  during  past  nuclear  tests  confirm  the  existence  of  substantial  ion- 

X  ^ 

ospheric  irregularity  following  high-altitude  detonations.  These 
irregularities  can  be  expected  to  reduce  system  performance  through  loss 
of  signal  coherence  across  both  frequency  and  time.  The  data  do  not,  how¬ 
ever,  suffice  in  themselves  to  provide  a  complete  characterization  of 
scatter  effects  at  HF  in  the  nuclear  environment.  Thus,  a  need  exists  for 
analytical  models,  based  on  a  general  understanding  of  such  effects,  that 
can  be  used  to  predict  these  effects.  This  final  report  describes  the 
work  performed  to  this  end  under  Tasks  1  and  2  of  Contract  DNA001-82-C-0218. 
Development  of  a  plan  for  a  field  experiment  to  complement  this  effort,  as 

called  for  by  Tasks  3  and  4  of  this  contract,  is  described  in  a  separate 

2 

report . 

Substantial  progress  has  been  achieved  in  recent  years  in  describing 
the  effects  of  ionospheric  irregularities  on  higher  frequency  signals. 


Superscripts  denote  references  listed  at  the  end  of  the  report. 
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such  as  those  used  in  transionospher ic  ground-to-satellite  communications 
links  and  by  missile-tracking  radar  systems.  As  a  result  of  this  progress, 
the  characteristics  of  the  irregularities  themselves  and  of  their  effects 
on  these  higher  frequency  signals  are  now  much  better  understood  for  both 
natural  and  nuclear  environments.  The  calculation  of  irregularity  effects 
is,  however,  inherently  simpler  at  the  higher  frequencies  than  it  is  at 
HF.  (Compare,  for  example,  the  analysis  in  References  3  and  4.)  Conse¬ 
quently,  the  modeling  of  scatter  effects  at  HF  has  received  only  sporadic 
attention  in  the  absence  of  new  insights  on  methods  for  their  calculation. 


Current  HF  nuclear-effects  propagation  codes,  such  as  NUCOM,  describe 
well  the  gross  multipath  effects  associated  with  ionospheric  refraction  at 
HF,  but  they  do  not  address  the  effects  of  scatter  from  ionospheric  irregu¬ 
larities.  Such  scatter  can  spread  the  signal  substantially  both  in  time 
and  frequency.  Our  approach  to  providing  a  complete  time-delay  and  Doppler 
spread  channel  characterization  at  HF  is  to  deal  separate!}  with  the  scatter 
effects  on  each  of  the  possibly  several  ray  paths  that  constitute  the  gross 
multipath.  Apart  from  the  strong  bending  of  the  HF  paths  by  refraction, 
each  of  these  paths  then  resembles  a  1 ine-of-sight  path  such  as  is  encoun¬ 
tered  in  transionospher ic  propagation  at  higher  frequencies. 

Recent  work  by  Flatte  et  al.^  on  the  scatter  of  oceanic  acoustic 
waves  has  shown  this  approach  to  the  determination  of  scatter  effects  in 
the  presence  of  strong  refraction  to  be  very  fruitful.  Furthermore,  as 
has  been  poi  ted  out  in  recent  JASON  studies^,  the  problem  addressed  by 
Flatte  et  al.  is  very  similar  to  that  of  the  scatter  of  obliquely  propa¬ 
gating  HF  radio  waves  by  ionospheric  irregularities.  The  ocean  and  the 
ionosphere  are  sufficiently  similar  media,  in  the  context  of  the 
refractive-index  variations  they  present  to  signals  propagating  obliquely 
through  them,  for  many  of  the  approximations  used  in  the  acoustic  analysis 
to  remain  valid  in  the  ionospheric  context.  Consequently,  adaption  of 
this  work  to  HF  radio  is  relatively  straightforward  and  yields  many  use¬ 
ful  results  without  inordinate  effort. 

The  results  of  the  analysis  undertaken  by  Flatte  et  al.^  also  bear, 


as  they  emphasize,  a  clear  relationship  to  the  scatter  models  that  have 
been  developed  for  1 ine-of -sight  propagation  through  slowly-varying  media. 


The  scatter  model  developed  recently  by  DNA'  for  transionospher ic  radio¬ 
wave  propagation  is  an  example  of  these  models.  Adaption  of  the  oceanic- 
scatter  theory  to  HF  thus  additionally  offers  the  prospect  of  a  unified 
approach  to  the  calculation  of  scatter  effects  in  the  nuclear  environment 
for  both  HF  and  higher  frequencies. 

A  preliminary  examination  of  the  scatter  theory  developed  by  Flatte 
et  al.^,  although  encouraging  regarding  the  above  expectations,  left  us 
uncertain  on  certain  points  as  to  the  approximations  involved  in  its 
development.  The  largely  hueristic  approach  adopted  by  Flatte  et  al. 
in  their  exposition  of  the  theory  relies  heavily  on  the  use  of  ray  con¬ 
cepts  to  extend  results  obtained  in  the  simpler  case  of  propagation 
through  statistically  uniform  media.  These  ray  concepts,  although  useful 
as  an  expository  tool,  become  suspect  in  the  very  regions--namely ,  near 
the  turning  points  of  rays  refracted  by  the  large-scale  nonrandom  varia¬ 
tions  in  the  background  medium--that  most  strongly  differentiate  the 
problem  from  that  of  a  uniform  medium.  Consequently,  we  felt  it  necessary 
in  this  initial  work  to  dwell  at  some  length  on  the  precise  manner  in 
which  these  regions  contribute  to  the  scatter  effects. 

This  question  is  examined  in  Section  3,  after  a  brief  introduction 
to  the  theory  developed  by  Flatte  et  al.',  as  adapted  by  us  to  HF,  in 
Section  2.  In  our  examination  we  employ  a  simple  ionospheric  model  that 
lends  itself  to  analytic  treatment.  Difficulties  associated  with  caustics 
formed  by  the  unscattered  ray  paths  are  analyzed,  and  an  improved  treat¬ 
ment  of  the  effects  of  these  caustics  is  described.  The  practical  impact 
of  these  effects  is  found  from  this  examination  to  be  slight. 

A  further  simplification  of  the  theory,  the  apex  approximation,  is 
discussed  in  Section  4.  This  approximation  attempts  to  exploit  the  ten¬ 
dency  of  scatter  to  occur  most  stronglv  near  the  apex  of  a  ray  path. 

The  e.xample  examined  bv  us  leads  us  to  conclude  that  this  approximation 
must  be  used  cautiously  at  HF. 

We  finally  consider  briefly,  in  Section  5,  the  r  ionship  of  the 
theory  to  the  channel  parameters  that  will  be  measured  experimentally. 
Expressions  are  given,  based  on  the  oce  n-scatter  results  of  Flatte  et  al. 


for  calculation  of  the  basic  parameters  that  an  experiment  would  be 
expected  to  measure.  These  results  are  presented  in  terms  of  the  two 
theoretical  quantities,  the  scatter  strength  and  diffraction  parameters, 
that  characterize  the  scatter  effects.  Although  the  results  presented 
here  suffice  to  demonstrate  the  completeness  and  relevance  of  the  theory, 
further  work  to  develop  expressions  tailored  to  HF  would  be  useful. 

Section  6  summarizes  our  results  and  suggests  an  approach  toward 
implementation  of  the  theory  for  HF  application. 


SECTION  2 

PATH- INTEGRAL  THEORY 


In  recent  work  on  the  scatter  of  acoustic  signals  by  oceanic  irregu¬ 
larities,  Flatte  et  al.^  have,  under  the  assumption  that  small-angle 
scattering  is  dominant,  expressed  the  effects  of  scatter  from  fluctuations 
in  the  refractive  index  of  the  propagation  medium  in  terms  of  an  integral 
over  all  possible  infinitesimally  segmented  ray  paths  connecting  the 
transmitter  and  the  receiver.  They  then  demonstrate  that  the  consequent 
"path- integral"  expressions  for  the  received-signal  characteristics  can 
be  reduced  to  approximate  expressions  involving  integrals  of  various 
quantities  taken  over  the  mean  (essentially,  unscattered)  ray  path  that 
results  from  the  large-scale  (and  deterministic)  spatial  variations  in 
refractive  index.  In  this  form,  the  expressions  that  describe  the 
scattered  field  bear  a  close  relationship  to  those  that  result  in  the 
simpler  problem  of  scatter  from  fluctuations  in  an  otherwise  homogeneous 
propagation  medium.  Extension  of  these  results  to  the  more  complex 
case  of  initially  refracted  unperturbed  ray  paths  entails,  as  will  be 
seen  below,  little  more  than  appropriate  generalization  of  the  definitions 
of  key  scatter  parameters.  Further  simplifications  that  can  be  made  in 
certain  cases,  such  as  the  thin-phase-screen  approximation,  also  have 
easily  recognized  analogs  in  the  theory  developed  by  Flatte  and  his 
col  leagues. 

As  is  discussed  in  detail  by  Flatte  et  al.,^  the  characteristics  of 
the  received  signal  are  determined  essentially  by  the  strength  and  the 
size  of  the  random  inhomogeneities  in  the  propagation  medium.  These 
properties  of  the  medium  are,  in  turn,  described  by  two  parameters.  The 
parameter  characterizes  the  strength  of  the  irregularities  in  terms  of 
the  correlation,  p(r,  r'),  between  the  phase-speed  fluctuations  at 
different  points  on  the  unperturbed  ray  path: 


(1) 


^  <  y.s  J 


ds  /  ds '  p( s , s  ■) 
ray  ray 


where  is  the  free-space  (electromagnetic)  wave  number  and  ds  is  an 
increment  of  distance  along  the  ray.  This  expression  constitutes  a 
straightforward  generalization  of  the  expression  that  arises  for  a  uni¬ 
form  background  medium. 

If  the  medium  fluctuations  decorrelate  over  a  distance  that  is  short 

compared  with  that  over  which  the  statistical  inhomogeneity  of  the  medium 

becomes  evident,  then  <t  can  be  written  in  terms  of  the  local  mean-squars 
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refractive-index  fluctuations,  <p  (s)>,  and  a  local  integral  scale,  L  (s) 
that  characterizes  the  spatial  correlation  of  the  fluctuations  in  the 
vi c ini ty  of  s : 
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For  a  horizontally  stratified  medium,  with  z  vertical  and 

ular  to  Che  planes  of  stratification,  the  dependence  on  s 

a  dependence  on  z,  and  L  can  be  written 
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In  Eq.  (  3 1 ,  u  is  measured  from  s  along  the  tangent  line  to  the  ray  at  s, 

at  angle  0  to  the  vertical,  as  is  Illustrated  in  Figure  1.  The  presence 

of  9  and  (j)  on  the  right  side  of  Eq.  (3)  indicates  that  Lp  depends  on  the 

direction  of  the  path  at  s  as  well  as  on  the  local  characteristics  of  the 

medium  there.  The  dependence  of  L  on  the  azimuthal  angle  ©  is  aosent  or 

^  ,  5 

the  oceanic  fluctuations  treated  bv  Flacte  et  al. 


where  n(k)  is  the  irregularity  wavenumber  spectrum 

fi(k)  =  (6) 

k" 

with  n  the  three  dimensional  spectral  index. 

The  parameter  A(s)  that  appears  in  Eq.  (4)  is  termed  the  phase  cur¬ 
vature.  It  is  defined  ’  as  the  second  derivative  of  the  phase 

path  length  for  a  lateral  displacement  of  the  ray  path  from  its  nominal 
position.  Thus,  A  is  related  to  the  Fresnel-zone  radius,  R  ,  by 
Rp.  =  2n(k^A)  .  Note  that  A  is  a  function  of  location  along  the  ray 
path,  as  is  Rp.  In  the  case  of  a  horizontally  stratified  medium,  the 
phase  curvature  for  horizontal  displacements  transverse  to  the  path  is 
the  same  as  that  for  a  straight-line  path  of  the  same  length  through  a 
uniform  medium  .  For  vertical  displacements,  however,  the  refractive 
bending  of  the  ray  path  markedly  affects  the  phase  curvature.  Addition¬ 
ally,  consideration  of  vertical  displacements  in  the  vicinity  of  caustics 
of  the  ray  path  raises  questions  that  involve  the  basic  limitations  of 
ray  theory  itself.  The  nature  of  these  difficulties  and  their  impact  on 
calculation  of  the  diffraction  parameter  is  discussed  in  Section  3. 

Evaluation  of  correlations  between  samples  of  the  received  signal 
additionally  requires  knowledge  of  the  phase-structure  function,  D, 
which  describes  the  difference  in  phase  between  two  points  separated 
in  space,  time,  or  frequency.  In  analogy  with  the  definition  of  D  for 
an  homogeneous  unperturbed  medium,  Flatte  et  al.^  define  this  quantity 
as 


D  = 


As  in  the  case  of  the  strength  parameter  '5,  the  generalization  of  the 
expression  for  D  to  a  refractive  background  medium  is  nominally  straight¬ 
forward  . 
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SECTION  3 

EFFECTS  OF  CAUSTICS 


In  Section  2  we  noted  that  the  complexity  of  the  unscattered  field 
structure  in  the  ionosphere  causes  some  difficulties  in  the  calculation 
of  the  diffraction  parameter  A.  These  difficulties  result  specifically 
from  the  presence  of  caustics  in  the  unperturbed  ray  structure.  Calcu¬ 
lation  of  the  diffraction  parameter  entails  determination  of  the  wave- 
field  characteristics  in  the  vicinity  of  the  unscattered  ray  path. 

Flatte  et  al.^  use  ray  expressions  to  describe  these  characteristics; 
their  expressions  fail  when  the  unperturbed  ray  is  near  a  caustic. 

The  ramifications  of  this  complication  can  best  be  assessed  in  the 
context  of  an  Ionospheric  model  that  has  analytic  expressions  describing 
the  unperturbed  wave  fields.  A  linear-gradient  model  of  the  variation 
of  ionospheric  electron  density  with  height,  although  somewhat  too  ideal¬ 
ized  for  practical  application,  serves  this  purpose  well.  This  model 
will  be  described  next,  after  which  the  nature  of  the  fields  near 
caustics  will  be  examined.  Modifications  to  the  calculation  of  the 
diffraction  parameter  that  account  for  the  caustic  effects  will  be 
described  in  the  final  portion  of  this  section. 

3 . 1  Linear-Gradient  Profile  Model 

A  linear-gradient  model  of  ionospheric  electron-density  variation 
with  height  allows  the  field  characteristics  to  be  determined  analytically 
With  this  model  the  variation  of  the  wave  fields  as  a  function  of  loca¬ 
tion  can  be  examined  to  whatever  degree  of  detail  is  required.  Basic 
expressions  for  the  ray  paths  and  their  phase-path  lengths  have  been 

g 

derived  from  the  general  formulas  given  by  Budden.  With  these,  the 
phase-path  length  (i.e.,  the  distance  along  a  ray  measured  in  wave¬ 
lengths)  to  a  scattering  point  can  readily  be  found  as  a  function  of  the 
position  of  this  point.  The  full-wave  characteristics  of  the  field  var¬ 
iations  near  caustics  can  also  be  determined  relatively  easily. 


The  linear-gradient  electron-density  profile  takes  the  form 


0  z  <  h 

oi(z  -  h  )  h  <  z 
o  o 


where  f^  is  the  plasma  frequency,  z  the  height,  h^  the  base  height  of 
the  profile  below  which  the  medium  is  free  space,  and  a  the  profile  gra¬ 
dient  parameter.  Ray  paths  are  described  generally  by 


o  r  <iz' 

X  =  S  /  - 

J  ^ 


where  x  is  the  horizontal  distance,  and  S  =  sin  4^  is  the  sine  of  the 

angle  between  the  wave-normal  (and  ray)  direction  and  the  vertical 

axis  at  the  origin  (x  =  0,  z  =  0).  The  quantity  q  is  the  refractive- 

9 

index  parameter  Introduced  by  Booker. 


q  =  n  COS0  =  (n^  -  )^  ,  (10) 

2  Z  h 

where  n=(l-f/f)is  the  refractive  index  for  frequency  f.  The 
P 

earth's  magnetic  field  has  not  been  included  in  the  refractive  index 
expression  because  its  effects  are  Incidental  to  those  of  primary  inter¬ 
est  here. 

The  phase-path  length  along  a  ray  from  the  origin  (0,0)  to  the 
point  (x,z)  is  given  generally  by 


Mi' 


<p  =  k(xS  +  /  qdz') 


where  k  =  2«f/c  is  the  free-space  wave  number.  The  point  (x,z)  may  be 
accessible  along  more  than  one  ray.  The  phase-path  length  ip  generally 
differs  in  this  case  for  the  different  rays,  as  is  implied  by  the 
presence  in  Eq.  (11)  of  the  parameter  S,  by  which  the  rays  are  di^'er- 
entiated  from  one  another. 


According  to  Eq.  (8),  the  ray  paths  within  the  ionosphere  are  sec¬ 
tions  of  parabolas;  beneath  the  ionosphere,  of  course,  they  are  sections 
of  straight  lines.  Any  point  beneath  the  ionosphere  is  connected  to  the 
source  by  two  rays.  One  of  these  is  the  straight,  line-of-sight  ray 
beneath  the  ionosphere;  the  other  ray  is  that  reflected  from  the  iono¬ 
sphere  by  the  process  of  refraction.  There  is  a  skip  zone  for  frequen¬ 
cies  sufficiently  high  to  penetrate  the  ionosphere  at  steep  angles  of 
incidence.  This  phenomenon  is  absent  for  the  linear-gradient  ionospheric 
model  because  the  electron  density  in  this  model  increases  without  limit 
with  increasing  height  in  the  ionosphere. 

Within  the  ionosphere,  a  shadow  boundary,  or  ray  caustic,  exists 
above  which  no  rays  penetrate.  The  location  of  this  boundary,  which  is 
evident  in  Figure  2,  is  frequency  dependent.  At  any  point  below  the 
shadow  boundary,  as  beneath  the  ionosphere,  at  least  two  rays  can  gener¬ 
ally  be  found  to  connect  a  given  point  with  the  source.  The  existence 
of  at  least  two  rays  can  be  deduced  qualitatively  by  considering  the 
evolution  of  the  line  of  sight  and  the  ionospherically  refracted  rays  as 
a  receiver  is  moved  up  into  the  ionosphere  from  beneath  it.  The  two  rays 
ultimately  coalesce  at  the  caustic.  Near  the  caustic,  as  can  be  seen  in 
Figure  2,  one  of  the  rays,  the  one  with  the  steeper  launch  angle,  has 
already  touched  the  caustic  at  a  slightly  shorter  range,  while  the  other 
ray  is  about  to  reach  the  caustic  at  a  slightly  greater  range. 

This  circumstance  has  implications  for  scatter;  multiple  ray  paths 
must  be  considered  between  the  scatter  location  and  both  the  source  and 
the  receiver  locations.  These  paths  are  illustrated  schematically  in 
Figure  3  for  a  scatter  point  at  an  arbitrary  point  on  unscattered  ray  path. 
If  the  two  possible  ray  paths  between  the  scatter  point  and  the  source,  or 
the  receiver,  constitute  distinct  rays  (a  condition  that  can  be  given 
precise  meaning),  Chen  only  one  of  them,  the  one  directly  evolvable  from 
the  unperturbed  ray  path  between  source  and  receiver,  is  important.  Near 
a  caustic,  however,  the  two  rays  are  not  distinct;  indeed,  a  ray  descrip¬ 
tion  of  the  field  behavior  near  the  caustic  fails  for  just  this  reason. 

5 

The  theory  developed  by  Flatte  et  al.,  which  uses  ray  expressions 
to  describe  the  fields  everywhere,  fails  to  address  the  above  situation. 
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Figure  3.  Multiple  rays  joining  point  to  path  endpoints 


The  hope,  of  course,  is  that  a  precise  treatment  of  shadow  boundaries 
formed  by  caustics  will  not  prove  to  be  critical.  This  expectation  is 
reasonable,  given  the  general  success  of  Kirchoff  approximations  in  the 
treatment  of  large  diffracting  objects,  but  such  a  result  should  not  be 
presumed  without  inquiry. 

Determination  of  the  locations  of  the  apex  and  caustic  surfaces  for 
an  ensemble  of  rays  facilitates  iteratively  selecting  the  appropriate 
ray  to  reach  a  particular  point.  This  procedure  will  be  used  to  examine 
the  rate  of  variation  of  some  path  parameters  with  displacement  of  a 
scattering  point  transverse  to  the  ray  path. 

For  our  linear-gradient  model,  the  apex  of  a  ray  represents  an 
extremum  of  the  height,  z,  as  a  function  of  the  distance,  x,  as  the  ray 
is  traced.  This  extremum  is  simply  identified  in  the  equations  for  the 
ray  paths,  [Eq.  (12)].  A  ray,  according  to  Eq.  (12),  reaches  its  apex 
at  x  =  D/2,  where  the  radical  in  the  expressions  valid  for  z  >  h  becomes 


Figure  4.  Locus  of  apexes  for  rays  traced  through  linear-gradient 
electron-density  profile 


Determination  of  the  location  of  the  caustic  surface  is  somewhat 
more  complicated.  A  caustic  surface  can  be  defined  as  the  locus  of 
points  2^)  for  which  the  rate  at  which  a  ray  is  shifted  laterally 

(transverse  to  itself)  by  a  change  in  its  launch  angle,  0^,,  is  zero. 
The  location  of  this  surface  can  most  simply  be  determined  for  the 
linear-gradient  profile  through  the  setting  to  zero  of  the  partial 
derivative  3x/00^  with  z  held  fixed,  with  the  result 


3.2  Field  Structure  near  Caustics 


Expressions  were  given  in  Section  3.1  for  finding  ray  paths  and  the 
phase-path  length  along  these  paths  from  a  localized  source  to  an  arbi¬ 
trary  point.  These  expressions  can  be  obtained  by  stationary-phase 
evaluation  of  Fourier-transf orm  integrals  that  relate  the  individual 

plane-wave  (angular  spectrum)  components  of  the  field  to  its  spatial 
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variation.  For  a  plane-stratified  background  medium  like  the  linear- 
gradient  ionospheric  profile  employed  here,  this  procedure  is  quite 
effective  and  yields  as  well  the  standard  ray-optics  expressions  for 
the  field  amplitude.  As  already  noted,  these  expressions  are  not  valid 
in  the  vicinity  of  caustics.  In  the  context  of  an  angular-spectral 
field  representation,  each  ray  is  associated  with  a  distinct  saddle 
point  across  which  the  path  of  the  transform  integration  passes.  The 
coalescence  of  rays  that  forms  a  caustic  is  thus  expressed  as  a  merging 
of  two  saddle  points,  and  the  stationary-phase  technique  for  approximat¬ 
ing  rhe  integral  across  the  saddle  as  an  asymptotic  series  fails.  The 
variation  of  field  intensity  near  caustics  was  first  determined  by  Airy; 
the  characteristics  of  this  variation  are  discussed  in  sufficient  detail 
for  our  needs  by  Budden  in  Section  11.7. 

The  variation  of  the  field-intensity,  I,  transverse  to  the  caustic 
surface  is  given  by 


I  =  Ai 


where  Ai  is  the  Airy  integral  function,  u  is  the  perpendicular  distance 
from  the  caustic,  k  is  the  local  wavenumber  in  the  medium,  and  R  is  the 
radius  of  curvature  of  the  caustic.  Equation  (16)  can  be  expressed  in 
terms  of  the  local  refractive  index  and  the  launch  angle  of  the  ray  that 
touches  the  caustic  at  u  =  o  through  use  of  Eq.  (10).  The  curvature  R 
can  also  be  expressed  in  terms  of  these  quantities  and  the  rate  at  which 
this  ray  moves  as  the  launch  angle  is  changed  by  noting  that  ''  is  given 


BigJ 


where  9  is  the  local  angle  of  the  ray  at  u=0.  Expression  of  R  in  terms 
o 

of  a  derivative  of  x  with  respect  to  the  angle  9  is  convenient  in  light 
of  the  explicit  dependence  of  x  on  9  given  for  the  linear-gradient  pro¬ 
file  by  Eq.  (12). 

Use  of  Eq.  (17)  in  Eq.  (16)  with  Eq.  (10)  employed  to  express  the 
local  wave  number  k  in  terms  of  the  free-space  wave  number  k^  through 
k  =  nk  ,  yields 

2  2 
2k  cos  9- 
o  I 

o 

2  .  2^  \3/2.2  |.-2, 

n  -sin  9j  j  d  x|d9 

o’  II 

o 

The  form  of  the  field  behavior  near  the  caustic  is  shown  in  Figure 
6,  in  which  the  behavior  of  the  Airy  function  is  shown  as  a  function  of 


Ai(x) 


Figure  6.  The  airy  integral  function  Ai(x) 


the  value  of  its  argument.  For  negative  values  of  the  argument,  the 

function  is  oscillatory.  If  the  Airy  function  is  expressed  as  an  Integra 

the  integral  can  be  well  approximated  by  stationary-phase  evaluation  at 

the  saddle  points  on  the  integration  path  for  sufficiently  large  values 

of  the  Airy-f unction  argument.  For  large  negative  values  of  the  argument 

this  path  crosses  two  saddle  points,  and  the  oscillatory  behavior  of  the 

function  can  be  interpreted  as  showing  interference  between  two  rays,  as 

their  relative  phase  changes,  on  the  illuminated  side  of  the  caustic. 

For  large  positive  values  of  its  argument,  Che  Airy  function  decays 
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In  Section  3.1  we  showed  that  the  ensemble  of  ionospheric  ray  paths 
forms  caustic  surfaces.  Ray  concepts  fail  to  describe  the  field  behavior 
near  these  caustics.  Thus  the  presence  of  caustics  in  regions  of  appre¬ 
ciable  scatter  can  greatly  complicate  attempts  to  calculate  scatter 
effects  using  a  purely  ray  treatment.  In  particular,  the  phase  curvature 
for  displacements  in  the  plane  of  the  unperturbed  ray  path  is  discontin¬ 
uous  at  caustics,  as  will  be  demonstrated.  The  effects  of  this  discon¬ 
tinuity  on  the  integral  for  the  diffraction  parameter  [Eq.  (4)]  is 
assessed  in  this  section. 

The  phase  curvature  is  defined^  as  the  second  derivative  of  the 
phase-path  length  with  respect  to  lateral  displacements  of  the  path, 
such  as  that  shown  in  Figure  7.  In  a  horizontally  stratified,  isotropic 

DISPLACEMENT 


UNPERTURBED 

PATH 


Figure  7.  Lateral  displacement  of  ray  path  vertically 


medium  (e.g.,  to  first  approximation,  the  ionosphere),  horizontal  dis¬ 
placements  yield  results  similar  to  those  for  a  homogeneous  medium. 
However,  displacements  in  the  plane  of  the  ray  path  lead  to  more  complex 
results  because  they  involve  the  phenomena  responsible  fcr  the  ray  refric 
t  ion . 
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If  the  phase-path  length  is  calculated  for  vertical  displacement  of 
the  ray  path  as  a  function  of  the  extent  of  the  displacement,  the  result 
varies  as  the  location  changes  on  the  path  from  which  the  displacement 
is  made.  The  results  of  several  such  calculations  are  shown  in  Figure 
8  for  vertical  displacements  at  a  series  of  locations  along  the  ray  path 
that  are  offset  horizontally  from  the  apex  by  multiples  of  2  km,  starting 
at  0  km  [Figure  8(a)].  The  phase  curvature  for  a  given  vertical  displace¬ 
ment  is  the  inverse  of  the  radius  of  curvature  of  the  curves  plotted  in 
Figure  8.  The  unperturbed  ray  path  in  these  examples  is  one  of  the  ray 
set  shown  in  Figure  2  for  a  linear-gradient  electron-density  profile  and 
a  5-MHz  frequency.  The  takeoff  angle  of  the  unperturbed  ray  was  adjusted 
to  a  value  that  gave  a  total  path  length  of  1500  km. 

The  phase-path  lengths  shown  in  Figure  8  for  the  perturbed  ray  paths 
were  calculated  as  the  sum  of  those  for  the  two  ray  segments  that  join 
the  displaced  path  point  to  the  two  (fixed)  endpoints  of  the  path.  The 
existence  of  multiple  phase-path  lengths  for  a  given  displacement  (i.e., 
the  existence  of  several  curves  in  each  of  the  figures)  results  from  the 
presence  of  more  than  one  ray  path  between  the  displaced  path  point  and 
one  or  both  of  the  path  endpoints.  Away  from  caustics,  only  one  pair  of 
these  paths  evolves  into  the  unperturbed  ray  path  as  the  vertical  di splace 
ment  is  reduced  to  zero.  This  unperturbed  ray  is  found  at  the  point  of 
zero  slope  in  the  variation  of  phase-path  length  with  vertical  displace¬ 
ment.  The  phase  curvature  used  to  determine  the  diffraction  parameter  A 
is  given  by  the  radius  of  curvature  at  this  point. 

At  the  ray  apex.  Figure  8(a),  the  unperturbed  ray  is  found  on  the 
lowest  of  the  three  curves,  where  the  phase-path  length  attains  a  maxi¬ 
mum.  The  other  two  curves  result  from  a  second,  higher  ray  on  one 
(middle  curve)  or  both  (upper  curve)  path  segments,  as  is  illustrated 
schematically  in  Figure  3  for  a  point  somewhat  to  one  side  of  the  ray 
apex.  (The  presence  of  multiple  rays  to  points  on  a  ray  path  can  also 
be  seen  in  the  ray  set  show^  in  Figure  2.)  As  the  path  is  displaced 
upward,  the  separation  between  the  iiiglier  and  the  lower  rays  decreases 
until  the  rays  coalesce  into  a  single  rav  at  a  caustic.  The  approach  to 
the  caustic  is  shown  in  Figure  8(a)  bv  the  decrease  in  the  difference  in 
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Figure  8.  Phase-path  lengths  as  a  function  of  bertical  displacement 
at  various  points  along  a  ray  path 


phase-path  length  among  the  three  curves  as  height  increases;  the  curves 
join  at  the  caustic  height. 

Figure  5,  introduced  previously,  shows  the  caustic  surface  pertinent 
to  the  ray  set  at  hand  here.  A  second  caustic  surface,  that  for  the  path 
with  the  s  :e  and  receiver  reversed,  also  is  pertinent.  This  surface 
can  be  found  by  mirroring  the  first  caustic  surface  in  a  vertical  line 
through  the  path  midpoint.  The  two  caustic  surfaces  thus  intersect  above 
the  path  apex. 

For  points  on  the  ray  path  horizontally  to  either  side  of  its  apex. 
Figures  8(b)  to  8(f),  the  middle  curve  of  Figure  8(a)  has  split  into  two 
curves.  This  split  results  from  the  asymmetry  between  the  two  sections 
of  the  path,  which  no  longer  traverse  equal  horizontal  distances.  For 
small  horizontal  offsets,  the  form  of  the  variation  in  phase  path  length 
as  the  vertical  displacement  changes  otherwise  is  generally  similar  to 
that  found  at  the  path  midpoint.  As  the  horizontal  offset  is  increased, 
however,  a  caustic  is  approached  on  the  ray  path,  and  the  junction  point 
of  the  lower  curve  pair  moves  nearer  the  position  of  the  undisplaced  ray 
on  the  lowest  curve  [Figures  8(d>  and  8(e)J.  As  the  horizontal  offset 
is  further  increased  so  that  the  caustic  is  passed  [Figure  8(f)]  the 
unperturbed  ray  shifts  to  the  upper  curve  of  this  pair,  with  an  attendant 

reversal  in  sign  of  the  phase  curvature.  This  sequence  is  shown  more 

clearly  in  Figure  9,  where  the  variation  of  phase  path  length  with  verti¬ 
cal  displacement  for  the  lower  curve  pairs  of  Figure  8  are  plotted  to  a 
slightly  different  scale  at  5-km  increments  in  horizontal  offset. 

The  variation  of  the  phase  curvature  along  the  unperturbed  ray  for 
a  path  of  about  1000-km  total  length  is  shown  in  Figure  10.  The  location 
on  the  ray  path  of  the  caustic  responsible  fo»-  the  discontinuity  evident 
in  the  curvature  at  about  450  km  is  shown  by  the  cross  mark  in  Figure  11. 

Only  one  half  of  the  ray  path  is  shown.  The  ray  traced  in  this  figure 

again  is  one  (that  launched  at  an  angle  of  75°  with  respect  to  the  zenith) 
of  the  set  shown  in  Figure  2. 

These  results  suggest  that  the  phase-path  length  might  generally  be 
considered  to  be  a  multivalued  function  of  the  position  of  a  point  joining 


Figure  10.  Phase  curvature  as  a  function  of  position  along  a  ray 


Figure  11.  Location  of  caustic  on  ray  with  75^  launch  angle 


rays  to  fixed  path  endpoints.  The  multivaluedness  represents  the  multi¬ 
plicity  of  rays  that  are  possible  for  each  of  the  path's  two  segments. 

At  caustics,  two  branches  of  the  phase-path  curve  join  as  two  rays  coalesce 
The  shift  of  the  phase  extremum  associated  with  the  total-path  ray  from 
one  to  the  other  of  these  two  branches  as  the  caustic  is  passed  causes 
the  phase  curvature  to  change  discontinuously ,  with  a  change  of  sign, 
across  the  caustic. 

This  multiplicity  of  ray  paths  is  of  little  consequence  for  scatter 
points  that  are  well  removed  from  caustics;  only  that  pair  of  path-segment 
rays  that  constitutes  a  valid  ray  for  the  entire  path  needs  to  be  con¬ 
sidered.  For  this  valid  ray,  the  phase  curvature  provides  a  useful 
estimate  of  the  Fresnel-zone  size,  which  in  turn  determines  the  distance 
to  either  side  of  the  ray  over  which  fluctuations  of  the  medium  must  be 
correlated  for  diffractive  effects  not  to  become  significant.  Other  com¬ 
binations  of  segment  rays  lack  the  extremum  of  phase  associated  with  a 
first-order  contribution  to  the  field  at  the  receiver. 

At  points  near  a  caustic,  however,  the  two  coalescing  rays  are  no 
longer  distinct  entities,  and  the  field  behavior  is  not  well  represented 
in  terms  of  rays  at  all.  The  Airy  integral  function  introduced  in  Section 
3.2  provides  a  better  description  of  the  fields  in  such  regions,  showing 
the  interference  between  the  two  rays  on  the  lighted  side  of  the  caustic, 
the  transition  into  the  shadow,  and  the  decay  of  the  field  in  the  shadow 
region.  Thus,  it  is  appropriate  to  seek  a  meaningful  length  scale  based 
on  this  description. 

Clearly,  refractive-index  fluctuations  well  within  the  shadow  region 
should  not  affect  the  received  field.  The  rate  at  which  the  field  decays 
as  the  shadow  region  is  entered  is  described  by  the  Airy  function;  this 
rate  provides  a  natural  length  scale  for  diffraction  near  the  caustic. 

This  characteristic  length  is  similar  to  the  distance  on  the  opposite 
(lighted)  side  of  the  caustic  required  for  the  difference  in  phase  be¬ 
tween  the  two  rays  to  become  equal  to  it  radians  (cf.  Figure  6).  Since 
the  ray  joining  the  perturbed  path  point  to  the  receiver  is  common  to 
both  rays,  this  latter  distance  is  also  that  over  which  the  difference 


in  total  phase-path  length  for  the  two  rays  becomes  n  radians  (as  seen, 
for  example,  in  Figure  9). 


Although  not  zero  at  the  caustic  (as  is  the  scale  derived  from  the 
phase  curvature),  the  characteristic  length  suggested  by  the  Airy-function 
behavior  remains  small  relative  to  the  size,  at  points  well  removed 
from  the  caustic,  of  a  Fresnel  zone.  The  contribution  to  the  diffraction 
integral,  Eq.  (4),  from  fluctuations  of  refractive  index  in  the  plane  of 
the  ray  path  near  the  caustic  consequently  would  appear  to  be  minimal, 
as  has  previously  been  assumed.^  Thus,  a  precise  treatment  of  the  caustic 
region  should  be  unnecessary,  and  use  of  the  phase  curvature  throughout 
the  diffraction  integral  would  not  be  expected  to  Introduce  substantial 
errors . 

The  Fresnel-zone  radius  given  by  the  phase  curvature  of  Figure  9 
is  shown  in  Figure  12  as  R^.  The  Fresnel-zone  radius  in  the  horizontal 
direction,  given  by  x(R-x)/k  R,  is  also  shown  for  comparison  as  R,  . 


Figure  12.  Fresnel-zone  radius  predicted  by  phase  curvature  as  a 
function  of  position  along  a  ray 

Between  the  caustics,  the  zone  size  in  the  vertical  direction  is  seen  tc 
remain  smaller  than  that  in  the  horizontal  direction.  The  diffraction 


SECTION  4 

APEX  APPROXIMATION 


The  effects  of  small-angle  forward  scatter,  resulting  from  stochastic 
irregularities  of  the  propagation  medium,  on  signals  that  also  are  refracted 
by  large-scale,  deterministic  variations  in  the  medium  can,  according  to 
Flatte  et  al.,^  be  described  to  a  good  approximation  in  terms  of  integrals 
of  various  quantities  taken  along  the  mean  (effectively,  unscattered)  ray 
path.  In  these  circumstances ,  the  intensity  of  the  smal ler-scale  random 
fluctuations  in  refractive  index  that  scatter  the  signal  may  be  related 
in  some  manner  to  the  deterministic  refractive-index  variations.  For 
example,  the  rms  magnitude  of  the  smaller  scale  electron-density  fluctua¬ 
tions  in  the  ionosphere  can  in  many  cases  be  approximated  as  being  some 
constant  fraction  of  the  mean  electron  density  at  each  point. 

A  further  approximation  to  the  path  integrals,  termed  the  apex 
approximation,  has  been  proposed  by  Flatte  et  al.^’  foj.  cases 

where  the  scatter  effects  are  concentrated,  as  in  the  above  example,  near 
the  apex  of  the  ray  path.  In  this  approximation,  the  random  phase  pertur¬ 
bations  accumulated  along  the  ray  path  are  applied  to  the  wavefront  in  a 
single  lumped  increment,  as  the  ray  passes  through  a  plane  placed  trans¬ 
verse  to  the  path  at  the  path  apex,  or  turning  point.  The  perturbed 
field  can  then  be  treated  as  a  secondary  source,  in  accordance  with 
Huygen's  principle,  to  calculate  the  received  wavefield  including  scatter. 
Thus,  the  apex  approximation  constitutes  a  generalization  to  refracted  ray 
paths  of  the  thin-phase -screen  approximation. 

5  S6C •  L 1  •  L 

Flatte  et  al.  ’  ’  "  note  that  the  validity  of  this  approxima¬ 

tion  requires,  in  addition  to  the  scatter  being  concentrated  near  the 
path  apex,  that  the  wave  intensity  in  the  vicinity  of  the  apex  not  be  much 
affected  by  the  scatter  and  that  the  phase  modulation  of  the  wave  be  cal¬ 
culable  using  geometrical  optics.  They  further  observe  that  these  require¬ 
ments  generally  are  met  if  propagation  remains  in  the  geometr ical -opt ic s 
regime  over  a  sufficient  horizontal  distance  about  the  apex,  namely 


1/2 

(r  *  ^  radius  of  curvature  at  the  apex  and  is 

the  vertical  scale  of  the  irregularities,  defined  in  Section  2.  This 
condition  implies  that  caustics  must  not  be  encountered  along  the  ray 
path  within  this  distance  since  geometrical  optics  is  fundamentally 
incapable  of  describing  the  wavefield  behavior  in  the  vicinity  of  a  caustic. 

Another  way  to  understand  this  constraint  is  to  consider  the  apex 
approximation  as  an  expansion  about  the  path  apex  of  the  parameters  that 
are  integrated  along  the  path  to  describe  the  scatter  effects.  In  order 
for  this  expansion  to  be  useful,  it  must  be  valid  for  the  entire  portion 
of  the  path  over  which  appreciable  scatter  occurs.  As  was  shown  in 
Section  3.3,  however,  the  phase  curvature  is  rather  ill  behaved  near 
caustics.  This  behavior  essentially  limits  the  range  about  the  path 
apex  over  which  the  phase  curvature  and  related  quantities  can  be  expanded. 

The  results  obtained  in  Section  3  provide  some  indication  of  the 
potential  usefulness  of  the  apex  approximation  at  HF.  Comparison  of  the 
caustic  surface  for  the  ray  set.  Figure  5,  with  the  apex  surface, 

Figure  4,  reveals  that  these  surfaces  are  not  greatly  separated  in  two 
regions,  namely  for  values  of  x  less  than  about  300  km  and  greater  Chan 
about  600  km,  corresponding  to  rays  with  very  steep  or  very  shallow 
take-off  angles,  respectively-  This  behavior  suggests  that  for  these 
rays  the  caustic  and  apex  locations  may  not  be  sufficiently  separated 
to  allow  the  use  of  the  apex  approximation. 

The  horizontal  distance  between  the  caustic  and  the  apex  can  be 
explicitly  calculated  using  Eqs.  (14a)  and  (15a)  to  determine  x^  and  x^ , 
respectively.  The  result  of  this  calculation  for  the  set  of  rays  presented 
in  Figure  2  is  shown  in  Figure  13.  Calculation  of  the  radius  of  curvature, 
r,  of  the  ray  path  for  a  linear-gradient  electron-density  profile  gives 


r  = 


2  2 

2f 


(19) 


at  the  path  apex.  Evaluation  of  Eq.  (19)  as  a  function  of  S  =  sin  0^, 
together  with  a  nominal  value  of  10  km  for  L^,  yields  the  second  curve 
shown  in  Figure  13.  Comparison  of  the  two  curves  confirms  that,  for  both 


caustic/apex  separation  with  apex-approximation  criterion 


the  very  steep  (0  <  25°)  and  the  very  shallow  (0^  >  75°)  rays,  the 

separation  of  the  path  apex  and  the  caustic  is  less  than  (r  x  L^) 
and  therefore,  according  to  the  criterion  reiterated  above.  Indeed  too 
small  for  the  apex  approximation  to  be  valid. 

Although  the  linear-gradient  profile  used  in  this  example  differs 
significantly  from  actual  ionospheric  profiles  in  some  of  its  character¬ 
istics,  the  shallow  rays  are  representative  of  the  propagation  behavior, 
on  long  HF  paths,  of  frequencies  below  the  maximum  usable  frequency.  Thus, 
it  is  necessary  to  conclude  that  the  utility  of  the  apex  approximation  may 
be  somewhat  limited  at  HF.  The  steep- incidence  rays,  although  representa¬ 
tive  of  the  behavior  on  a  short  path  of  frequencies  below  the  vertical- 
incidence  critical  frequency,  are  less  instructive;  they  represent  a  case 
for  which  the  basic  applicability  of  the  path-integral  scatter  theory  can 
be  questioned  more  generally.  Namely,  it  is  difficult  to  view  the  electron 
density  fluctuations  as  leading  only  to  small-angle  scatter  in  this  case 
since  the  signal  frequency  approaches  the  plasma  frequency  near  the  path 
apex,  and  the  fluctuations  consequently  strongly  perturb  the  height  at 
which  the  plasma  is  overdense  at  the  wave  frequency. 


SECTION  5 

CHANNEL-PROBE  PARAMETERS 


An  experiment  has  been  planned  to  measure  the  impulse  response 
(i.e.,  the  time-delay  and  Doppler  spread  of  the  received  signal  for  an 
impulsive  source)  of  an  HF  channel.  This  measurement  suffices  to  char¬ 
acterize  the  channel  completely  in  the  absence  of  spatial  diversity. 

The  theory  discussed  in  preceeding  sections  of  this  report  characterizes 
the  impulse  response  of  the  channel  in  terms  of  two  parameters,  the 
coherent  bandwidth  and  the  coherence  time,  that  are  closely  related  to 
the  parameters  that  are  measured  experimentally. 

The  coherent  bandwidth,  (n  )  is  given  by^’ 

w 

2  /  2  \-l 

^  u)  lAli^ln*)  (20) 


in  the  fully  saturated  case,  neglecting  a  factor  of  order  unity  that 
involves  the  ratio  of  a  scale  for  variation  of  the  phase-structure 
function  to  that  for  the  refractive-index  fluctuations.  The  delay  spread 
(Tp)  measured  experimentally  is  related  to  this  bandwidth  by 


T 


D 


Zx 
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The  coherence  time  characterizes  the  correlation  over  time  of  the 
signal  envelope  fluctuations.  The  latter  is  given  by^’ 


<'f  (t+At  )'l'(t  )>  =  exp  C-%D(At)] 


(22  ) 


where  f(t)  is  the  (complex)  signal  envelope  amplitude  (also 

reduced  wavefunction )  at  time  t  and  D  is  the  phase-structure 

defined  by  Eq.  (3).  Generally,  the  phase-structure  function 

5,  §7.5  .  ^  ^ 

time  separations  can  be  written  in  terms  of  the  rate 


called  the 
function 
for  small 
of  change 


of  the  strength  parameter  <t  as 


D(At)  =  (At)^*^  ,  (23) 

although  for  small  spectral-index  values  this  quadratic  approximation 
can  break  down. 

In  the  case  of  HF  ionospheric  propagation,  $  usually  is  considered 

to  be  determined  by  convection  of  irregularities  across  the  ray  path- 
6 

In  this  case, 

D(At)  ~  ,  (24) 


and  the  coherence  time,  T  ,  is 

c  ’ 

T 

c 


(25) 


where  v  is  the  convective  speed  of  ionospheric  winds  transverse  to  the 
path. 

Experimentally,  the  coherence  time  is  determined  by  measuring  the 
time  variation  of  a  CW  signal.  The  result  of  this  measurement  usually 
is  transformed  to  yield  a  Doppler  spectrum.  The  Doppler-spread  parameter 
(Dp  that  characterizes  the  width  of  this  spectrum  is  related  to  the  co¬ 
herence  time  by 


SECTION  6 

CONCLUSIONS  AND  RECOMMENDATIONS 


The  path- integral  theory  developed  by  Flatte  et  al.  to  describe 
scatter  effects  in  the  presence  of  deterministic  refraction  should  be  a 
useful  tool  in  HF  radiowave  applications.  A  major  simplification  is 
achieved  in  their  work  by  the  approximation  of  the  path- integral  expres¬ 
sions  by  integrals  taken  along  the  unscattered  ray  path.  We  have  examined 
carefully  the  effects  on  these  expressions  of  caustics  encountered  midway 
along  a  ray  path,  which  the  theory  as  presented  by  Flatte  et  al.  does  not 
address.  (Flatte  et  al.  do  note,  however,  that  their  results  are  not 
valid  if  the  receiver  is  located  near  a  caustic,  such  as  would  be  the 
case  for  an  HF  receiver  that  happened  to  be  near  the  skip  distance.)  Our 
examination  found  caustics  encountered  in  midpath  not  to  have  much  impact 
on  the  scatter  behavior.  Thus  the  inclusion  of  their  effects,  although 
fundamentally  correct,  would  represent  an  unnecessary  elaboration. 

The  apex  approximation,  which  reduces  the  calculation  of  scatter 
effects  on  a  strongly  refracted  ray  path  to  a  thin-phase-screen  model, 
constitutes  a  further  simplification  of  the  path-integral  theory.  This 
simplification  is  possible  only  if  caustics  encountered  along  the  path 
are  sufficiently  removed  from  the  path  apex.  Our  comparison,  for  a 
specific  example,  of  the  separation  between  caustic  and  apex  with  the 
minimum  distance  necessary  to  permit  use  of  the  apex  approximation  found 
this  condition  not  to  be  met  for  a  significant  range  of  ray  take-off 
angles.  Although  our  results  constitute  only  a  single  example,  they 
suggest  that  the  apex  approximation  should  be  used  cautiously  at  HF. 

As  was  noted  in  the  introduction  to  this  report,  the  linear-gradient 
ionospheric  model  used  in  this  work  falls  somewhat  short  of  the  degree  of 
realism  desired  for  practical  use.  A  better  ionospheric  model,  which 
still  retains  the  advantage  of  yielding  explicit  analytic  expressions 
for  the  unscattered  ray  path,  is  that  of  a  parabolic  variation  of  electron 
density  with  height.  This  model  probably  represents  the  simplest  that 


should  be  considered  for  practical  application.  We  recommend  that  the 

path-integral  scatter  theory,  as  reduced  by  Flatte  et  al.^  to  ray-path 

integrations,  be  implemented  using  this  ionospheric  model.  In  this 

implementation  expressions  for  the  basic  channel  parameters  should  be 

developed  specifically  for  the  HF  application.  Any  need  for  further 

elaboration  of  the  theory  (including  the  consideration  of  more  elaborate 

ionospheric  models)  can  better  be  addressed  after  results  obtained  from 

the  theory  in  this  form  have  been  compared  with  data  from  a  field  exper- 

2 

iment,  such  as  chat  described  in  the  companion  report  prepared  under 
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